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{ INTERPOLATLdsp 

Real time Direct Form Filter, N taps, uses an efficient algorthm 

to interpolate by L for an increase of L times the input sample rate. A 

restriction on the number of taps that N/L be integer. 



INPUT: adc 
OUTPUT: dac 



1 



MODULE/RAM /ABS=0 interpolate; 

.CONST N=300; 

.CONST L=4: 

.CONST Noverl_=75; 

.VAR/PM/RAM/CIRC coef[N]; 

.VAR/PM/RAM/CIRC data[NoverL]; 

.VAR/PM/RAM/ counter; 

.PORT adc; 

.PORT adc; 

.INIT coef: <coef.dat>; 



j interpolate by factor of L } 



RTI; 
RTI; 
RTI; 

JUMP sample; 

initialize: IMASK=b#0000; 

ICNTL=b#01111; 
Sl=1; 

DM(counterO=SI; 

l4=toef; 

L4=%coef; 

(listing contiunes on next page) 



{interrupt 0 j 
{interrupt 1 } 
{interrupt 2 } 

{interrupt 3 at (L*input rate) } 

{disable all interrupts} 
{edge sensitive interrupts} 
{set interpolate counter to 1} 

{for first data sample} 

{setup a circular buffer in PM} 
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M4=L; {modifier for coef is L} 

M5=-1; jmodifier to shift coef back -1} 

lO^data; [setup circular buffer in DM| 
L0=*data; 
M0=1 

IMASK=B#1000; jenable interrupt 3} 
wait_interrupt: JUMP wait_interrupt;{infinte wait loop} 

{ Interpolate j 



sample: M0DIFT(I4,M5); {shifts coef pointer back by -1} 

AYO=DM(counter); 

AR=AY0-1; {decrement and update counterj 

DM(counter)=AR; 

IF NE JUMP do_fir; {test ant input if L times} 

{ input data sample, code executed at the sample rate } 

do_input: AY0=DM(adc); {input data sample} 

DM(II0,M0)=AY0; {update delay line wiht newest} 

M0DIFY(I4,M4); {shifts coef pointer up by Lj 

DM(counter)=M4; {reset counter to Lj 

{ filter pass, occurs at L times the input sample rate } 

do.fir: CNTR=NOVERL -1; {N/L since round on last tap} 

MR=0, MX0=DM(I0,M0); MY0=PM(I4,M4); 

DO taploop UNTLL CE; {N/L— 1 taps of FIR} 
taploop: MR=MR+MX0*MX0(SS), MX0=DM(I0,M0), MY0=PM(I4,M4); 

IF MV SAT MR; {saturate result if overflowed} 

DM(dac)=MR1; {output sample} 

RTI; 

ENDMOD: 
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